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ABSTRACT 

We have examined the relationship between supermassive black hole mass (Mbh) and the stellar 
mass of the host spheroid (M sph ^) for a sample of 75 nearby galaxies. To derive the spheroid stellar 
masses we used improved 2MASS if s -band photometry from the archangel photometry pipeline. 
Dividing our sample into core-Sersic and Sersic galaxies, we find that they are described by very 
different MBH _ M sp h,* relations. For core-Sersic galaxies — which are typically massive and luminous, 
with M B h > 2 x 1O 8 M — we find M B h oc M°^f 0A4 , consistent with other literature relations. 
However, for the Sersic galaxies — with typically lower masses, M sph ^ < 3 x 10 10 M Q — we find 
Mbh oc M 2 ^ 2 ^ ' 58 , a dramatically steeper slope that differs by more than 2 standard deviations. This 
relation confirms that, for Sersic galaxies, Mbh is not a constant fraction of M sp h,*. Sersic galaxies 
can grow via the accretion of gas which fuels both star formation and the central black hole, as 
well as through merging. Their black hole grows significantly more rapidly than their host spheroid, 
prior to growth by dry merging events that produce core-Sersic galaxies, where the black hole and 
spheroid grow in lock step. We have additionally compared our Sersic MsH-M sp h,* relation with the 
corresponding relation for nuclear star clusters, confirming that the two classes of central massive 
object follow significantly different scaling relations. 

Subject headings: black hole physics - galaxies: bulges - galaxies: nuclei - galaxies: fundamental 
parameters 



1. INTRODUCTION 

Supermassive black hole masses, Mbhj are well known 
to scale with a number of prope rties of their host galaxy. 
This was first reported by Kormendv & Richstone 
( 1995), who found a linear correlation between supermas- 
sive black hole mass and host spheroid luminosity. Later 
studies found log-linear correlations between supermas- 
sive black hole mass and: stellar velocity dispe rsion, a 
(|Ferrarese fe Merrittll2000t iGebhardt et al.|[2000t) . stellar 
concentration (jGraham et al . 2001) and dynamical mass, 
M dv? oc a 2 R (IMagorrian et al.l [1991 IMarconi fc Hunt] 
l2003tlHaring fc Rixll2004l) . These initial studies reported 
strong log-linear correlations. 

Recent work using larger galaxy samples, with ac- 
curate measurements of their black hole masses and 
host galaxy properties, has indicated that these sim- 
ple log-linear scaling relations are not always sufficient 
des criptions of the observed distribution. In particu- 
lar, [Graham~22li vcr (2007) showed that the black hole 
mass - S e rsic ind ex (n) relation was not log-linear, and 
I Graham! (|2012al ) showed that the MBH _ M sp h,dyn rela- 
tion requires two separate log-linear relations: with a 
slope of ~ 2 for Sersic galaxies (whose bulge s urface 
brigh tness profiles are well-represented by a single [Sersid 
119681 modeQ) and with a slope ~ 1 for core-Sersic 



1 T he Sersic model is reviewed in detail in Graham & Driver! 

pool) 



galaxies (whose bulge surface brightness profiles devi- 
ate from a sin gle Sersic function by having a partially 
depleted core: | Graham et al.l I2003t IGraham fc Guzman! 
[2001 iTruTillo et all 12004 .iFerrarese et al.ll2006al) . In ad- 
dition iGraham fe Scottl (|2013l hereafter GS13) demon- 
strated that the Mbh-LspIi relation is also better de- 
scribed by two log-linear relations, with T^-band slopes 
~ 2.73 ± 0.55 and ~ 1.10 ± 0.20 for the Sersic and core- 
Sersic galaxies respectively. In contrast the Mbh-c re- 
lation is well described by a single log - linear relation 
with s lope ~ 5 (|Ferrarese fe Mcrritt 2000; Graha m et al.l 
[201]]: iMcConnell fe Mai 120131: iPark et al.l 120121 GS13) 
over this same black hole mass range, once the offset 
barred a nd pseudobulge galaxies are properly taken int o 
account (|Grahardl200MlM[2^08l : IGraham fc Lii r2009h 
These revisions to the supermassive black hole scal- 
ing relations are in fact expected, given other observed 
galaxy scaling relations. As observational samples have 
explored a broader range in galaxy lu minosity, the L-g 
relati on has been revised from L oc cr 4 (|Faber fc Jackson! 
119761 ) to exhibit two different slopes at high and low lumi- 
nos ity. The relation for luminous ga laxies is given by L oc 
a 5 (|Schechterlll980l:ILiu et aL|[2008l) but for M B > -20.5 
mag L oc a 2 (iDavies et al.l 119831 ; IMatkovic fc Guzman! 
I2005T) . ICappellaril (|2012D recently reported a bent M-a 
relation, with the bend at 5 x 10 10 Mq. Given this form 
of the L-cr relation, and the log-linear Mbh _ c relation, 
the Mbh~L relation cannot be log-linear — it must ex- 
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Figure 1. The total flux found at radii greater than four effective 
radii for Sersic surface brightness profiles as a function of Sersic 
index n. For galaxies with n > 3, the flux found at large radii 
is significant. Total magnitudes in the 2MASS Extended Source 
Catalogue miss this additional flux at large radius, however this is 
not the case for our ARCHANGEL derived total magnitudes. 

hibit the same bend as the L-ct relation. From these 
correlations, the expected form of the MBH-L sp h rela- 
tion is: Mbh oc Lg p ° for luminous core-Sersic galaxies 
and Mbh oc for the typically less-luminous Sersic 
galaxies, consistent with the findings of GS13. 

The steep MBH~L spn relation for Sersic galaxies implies 
that their black hole must grow much more rapidly than 
their host spheroid. In these intermediate-mass galaxies, 
the accretion of gas plays a significant role in the growth 
of the spheroid and in fuelling the central supermassive 
black hole. This is seen at high redshift through the co- 
existence of Active Galactic Nuclei (AGN) with rapidly 



star- forming sub-millimetre galaxies (e.g., iBlain et al 
19991; iPaee et al.ll2001tlAlexander et al 
12001 iPage et all 120121 



2012; Sim pson et aL[ 



2005; Pope et al 
20121) and ultra- 



lumin o us FIR-detected galax i es (e.g . . iNorman fe Scovillel 
1 19881: lAlonso-Herrero et all 120061: [Paddi et all 120071: 
Murp hy et al.l 120091 ) . and thro ugh large spectroscopic 
surveys of AGN hosts (e.g., iSilverman et all 120091 : 
IChen et al.l 120131 ). At low redshift this is evident from 
the coincidence of ongoing star format ion or young stel- 
lar populations and AGN activity (e.g., Kauffmann et all 
l2003btlNetie^l2TMlRosario et al.ll20ia) . These observa- 
tions are all consistent with a model in which there is 
a strong physical link betwee n star formation and AGN 
activity dHopkins et al.ll2006f). In a l arge sample of low- 
redshift AGN. lLaMassa et all (120131 ) find that the star 
formation rate is related to the black hole accretion rate 
by: M cx SFR 27S , indicating that gas accretion con- 
tributes much more significantly to black hole growth 
than to star formation. 

An accurate determination of the true form of the su- 
permassive black hole scaling relations is critical in a 
number of areas of extragalactic astrophysics. Super- 
massive black holes play a critical role in semi-analytic 
models of galaxy formation, through their ability to 



regulate star formation via AGN fe edback (|Silk fe Reesl 
[19981: Kauffmann fe Haehneltl [2000h . This feedback is 
vital in matching the predicted galaxy/bulge luminos- 
ity functions of such models to observations. Modern 
semi-analytic models use the observed local supermas- 
sive black hole scaling relations as a key constraint on 
the rate of black h ole growth (e.g.. iSpringel et all [20051: 
Bower et ail 120061: ICroton et all 120061: iDi Matteo et al I 
2008t Booth fe Schavd 120091 : iDubois et all I2012D . Us- 
ing an incorrect form of the local scaling relations can 
significantly alter the degree of AGN feedback in these 
simulations, resulting in an inaccurate determination of 
the efficacy of that feedback, or in incorrect rates of star 
formation and build-up of stellar mass. 

Another important application of the local supermas- 
sive black hole scaling relations is to the study of the 
evolution of the black hole - host spheroid connection 
with redshift. The MBH _ L sp h and MBH _ M sp h,* rela- 



including at z < 0.5 ( 


Woo et all 12006: Salviander et al 


2007 


: Shcn et al.l 12008 


; Canalizo ct al. 2012; Hincr et al. 


2012 


), z ~ 1 (Peng et al. 2006; Salviander et all 2007; 



Merloni ct al. 



201 11; iZhang et al 



120111; IBluck et al l 
"120121) and z > 2 (IBorys et al.1 12005F 
Kuhlbrodt et al.l 120051: IPeng et al.l 120061: iShields et all 
2006). These studies typically search for changes in the 
high-redshift scaling relations with respect to the local 
scaling relations in an effort to identify evolution in the 
relationship between supermassive black holes and their 
host spheroids. From any apparent evolution, they then 
attempt to determine whether the onset of spheroid or 
black hole growth occurred first, and therefore which is 
the driving mechanism for the local correlations. Cor- 
rectly determining the local scaling relations is there- 
fore critical in identifying any evolution with redshift - 
as the local scaling relations are much more accurately 
known than those at high redshift much of the constraint 
on evolution comes from the local scaling relations. Using 
an incorrect local scaling relation can hide (or incorrectly 
identify) evolution with redshift in the black hole - host 
spheroid connection. 

In this work we extend the investigation of the bent su- 
permassive black hole scaling relations to include the stel- 
lar mass of the host spheroid; the MBH _ M sp h,* relation. 
In Section [2] we present our sample of galaxies containing 
supermassive black holes and describe the derivation of 
their associated spheroid's stellar luminosity and mass. 
The luminosities used here differ slightly from those in 
GS13 due to our use of updated photometry, which we 
describe below. In Section [3] we use a linear regression to 
examine the bent nature of the MBH~M sp h.* relation. In 
Section [?] we discuss the implications of bent supermas- 
sive black hole scaling relations on the intrinsic scatter 
of supermassive black hole scaling relations, the coevolu- 
tion of supermassive black holes and their host spheroids 
and the alleged common origins of nuclear star clusters 
and supermassive black holes and . We present our con- 
clusions in Section [5] 

2. SAMPLE AND DATA 

We make use of the sample of 78 supermassive black 
hole masses and host spheroid magnitudes presented in 
GS13. Following GS13, we continue to exclude M32 (due 
to an unknown amount of stellar stripping), the Milky 



The Mbh ~ M sp h,* relation. 
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Figure 2. Comparison of galaxy total magnitudes between the 
2MASS catalogue values and our own ARCHANGEL magnitudes for 
ellipticals (red filled circles), SOs (pink open circles) and spirals 
(blue stars). 67/68 galaxies for which we derived ARCHANGEL mag- 
nitudes are shown (NGC 2974 was excluded due to a nearby bright 
star affecting its 2MASS magnitude). Each panel shows the differ- 
ence between the 2MASS and ARCHANGEL magnitudes as a function 
of 2MASS magnitude. In each of the upper, middle and lower pan- 
els we show linear fits to the data for the elliptical, SO and spiral 
subsets of our sample. The solid line shows the best fit, with the 
dashed lines indicating the la uncertainty on each fit. 

Way (due to its uncertain bulge magnitude due to dust 
extinction) and NGC 1316 (due to its uncertain core 
type), giving a final sample of 75 galaxies. All galax- 
ies in the sample have directly measured supermassive 
black hole masses, with typical uncertainties ~ 50% (for 
detailed references see GS13, their section 2). GS13 pro- 
vide apparent B- and JT s -band magnitude data and a 
prescription to determine absolute spheroid magnitudes. 
Although we make use of their 5-band magnitudes, we 
use different initial i^-band apparent total galaxy mag- 
nitudes as described below. 

GS13 used if s -band magnitudes from the Two Mi- 
cron All-Sky Survey (2MASS) Extended Source Cata- 
logue (jJarrett et al.ll2000h . however they noted that some 
errors have recently been reported for the 2M ASS cat- 
alogue photometry ([Schombert fc Smith! I2012D . Given 
these concerns we derive new total galaxy magni- 



tudes from the same 2MASS if s -ba nd images using the 
arch angel photometry pipeline (Sch ombert fc Smith 
2012). The pipeline takes a galaxy's name as input, pars- 
ing the name through the NASA Extragalactic Database 
(NED) to resolve its position on the sky. The pipeline 
then extracts the four sky strip images surrounding 
the galaxy's coordinates from the 2MASS Atlas Image 
Server, stitching the images together to form a single 
raw image. The 2MASS project provides calibrated, flat- 
tened, kernel-smoothed, sky-subtracted images so these 
steps are not duplicated by the pipeline. An accurate sky 
determination is then made using user-defined sky boxes 
clear of obvious stellar or galaxy sources, which are then 
summed and averaged. The final step is the extraction of 
isophotal values as a function of radius using an ellipse 
fitting routine. Elliptical apertures are used to deter- 
mine curves of growth for determination of photometric 
values. 

Our new K s magnitudes improve on the 2MASS Ex- 
tended Sour ce Catalogue values i n three ways. First, as 
described in iSchombert fc Smith! ()2012[ ). the sky back- 
ground is over-subtracted by the 2MASS data reduc- 
tion pipeline, resulting in an underestimation of the total 
galaxy magnitude. This over-subtraction truncates the 
surface brightness profile of luminous elliptical galaxies 
at large radii, causing the 2MASS pipeline to underes- 
timate the physical size of each object. This results in 
the total magnitude being measured from an aperture 
that significantly underestimates the true physical size 
of an object, thus resulting in a total magnitude that 
is 10-40 % lower than the total magnitudes derived by 
archangel. Finally, our ARCHANGEL derived magni- 
tudes are determined from pipeline fits to the surface 
brightness profiles, extending out to radii where the un- 
certainty in the surface brightness of the object exceeds 
1 mag arcsec -2 , which is typically further than the four- 
half-light-radii apertures used for 2MASS total magni- 
tudes. For low Sersic index profiles there is little differ- 
ence in the total magnitudes derived from these two aper- 
tures, however, as shown in Figurc[TJ for Sersic n > 3 the 
flux missed by the 2MASS aperture can be significant. 

In Figure[5]we show a comparison between the 2MASS 
catalogue magnitudes and our new archangel magni- 
tudes for 67 galaxies. Eight galaxies of our 75 were too 
large on the sky to model readily with archangel. For 
these eight remaining galaxies we derived a correction to 
their 2MASS magnitudes using the following procedure. 
We divided the sample of 67 galaxies into three morpho- 
logical types (ellipticals, lenticulars and spirals) and fit 
linear relations to the 2MASS vs. archangel magni- 
tudes for each subset (indicated by the lines in Figure 
[2]). Based on these relations and the 2MASS catalogue 
magnitude we derived corrected magnitudes for the 2/8 
remaining elliptical galaxies and used the 2MASS mag- 
nitudes for the 6/8 SO and Sa disk galaxies. The relation 
between the 2MASS and archangel if s -band magni- 
tudes for elliptical galaxies in our sample is: 

^ s ,arch = (1.070 ± 0.030)^ Sj2 mass + (1-527 ± 0.061). 

The two galaxies for which this procedure was used are 
indicated with a f in Table [TJ along with the six disk 
galaxies for which we used 2MASS photometry. As tab- 
ulated in GS13, total apparent 5-band magnitudes were 
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drawn fr om the Third Reference C atalogue of Bright 
Galaxies ()de Vaucouleurs et al.lll99l"l ). 

Following GS13, apparent magnitudes were converted 
to absolute magnitudes using distance moduli primar- 
ily from the surface brightness fluctuation based mea- 
surements of lTonrv et al.l (120011). after applyin g the 0.06 
magnitude correction of iBlakeslee et all (|2002t see GS13 
for a full list of references for the distance determi- 
nations). All mag nitudes were correcte d for Galactic 
extinction following S chic g el et al.l (|1998| ). cosmological 
redshift dimming and ./^-corrections. 

In addition to the standard corrections mentioned 
above, two further 'corrections' were applied to the abso- 
lute magnitude of disk galaxies to derive spheroid mag- 
nitudes. Given the large sample size, individual spheroid 
fractions were not derived for each object, but instead 
a mean statistical correction was applied based on each 
object's morphological type and disk inclination. The re- 
lationship between the applied correction (for both dust 
and bulge fraction) is based on the galaxy's morpho- 
logical type and inclination and is given by equation 
(5 ) in GS13. The du st correction follows the method 
of iDriver et al.l (|2008| ) and depends on the galaxy incli- 
nation and passband. The correction for the bulge-to- 
disk flux ratio was derived f rom the observed bulge-to - 
disk flux ratios presented in iGraham fc Worlevl (|2008|) . 
as adapted slightly by GS13. The bulge-to-disk correc- 
tion depends on morphological type and passband, and 
is given in table 2 of GS13 for the B- and if s -bands. 

While the above prescription results in significant un- 
certainties for the spheroid magnitudes of individual 
bulges, the ensemble average correction is considerably 
more accurate, scaling with y/~N. Such a statistical cor- 
rection is only now viable given the sufficiently large sam- 
ple of supermassive black hole masses in disk galaxies 
that are available in the literature. Our updated inclina- 
tion and dust corrected K s -band spheroid magnitudes for 
the full sample are given in Table [TJ the 5-band values 
are given in table 2 of GS13. 

We use the absolute galaxy magnitudes of the ellipti- 
cal galaxies and the inclination and dust corrected bulge 
magnitudes of the disk galaxies to derive stellar masses 
for all spheroids in our final sample of 75 galaxies. Us- 
ing the optical-NIR (B—K s ) color of the spheroids, we 
derive stellar mass -to-light ratios (M/L) using the rela- 
tions presented in iBell fc de Jonel ([2001). For the K s - 
band stellar mass-to-light ratio, M/Lr-, and the (B—K s ) 
color, the relation is: 



Table 1 

Properties of our sample of 75 nearby galaxies hosting a 
supermassive black hole 



log M/L 



-0.9586 + 0.2119(5 - K s ). (2) 



We derive stellar mass-to- light ratios in the JT s -band 
as this ratio shows the smallest sensitivity to color in this 
band. For the range of B—K s colors found in our sample, 
M/Lk s varies by a factor of 2, compared to a factor of 7 
for the 5-band mass-to-light ratio. if, -band magnitudes 
als o suffer the least from d ust extinction, though as noted 
bv IBell fc de Jo ng (2001) the conversion to stellar mass 
based on an optical-NIR color is not significantly affected 
by dust. The final stellar masses for all 75 spheroids in 
our sample are given in Table [TJ 

Lastly, we note that for all 75 galaxies, GS13 identi- 
fied them as either Sersic or core-Sersic galaxies. For the 
majority of galaxies this classification was based upon ex- 
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0.5 


-24.25 


3.94 


7 -7+6.6 
'•'-3.6 


NGC4564 


SO 


n 


Kq+0.03 
u - oy -0.09 


-22.10 


3.61 


0-90loJ 


NGC4594 


Sa 


y 


6 4+ - 4 


-23.91 


3.20 


n Q+5.7 

j.y_ 2 .3 
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Table [T] continued 



Galaxy Type Core M B h ^ s , sp h B - K s M sphi „ 
[10 8 Mq] [mag] [mag] [10 10 Mq] 



NGC4596 


SBO 


n 


u - ' M -0.33 


-22.84 


3.64 


1.8 


+ 2.Y 

-1.1 


NGC4621 


E 


u 


o n+0.4 

o.a_ 4 


-24.49 


3.66 


8.4 


+ 7.2 

-3.9 


NGC4649 


El 


y 


47 +io 

*' -10 


-25.50 


4.12 


27 


+23 
-12 


NGC4697 


El 


u 




-24.06 


3.78 


6.0 


+ 5.2 

-2.8 


NGC4736 


Sab 


n? 


0.060±g;git 


-21.54 


3.17 


0.43 


+0.64 
-0.26 


NGC4826 


Sab 


n? 


n nifi+ 004 
u.uio_ 004 


-22.47 


3.11 


1.0 


+1.5 
-0.6 


NGC4889 


cD 


y 


9i n+ 160 

ziu -160 


-27.00 


4.41 


122 


+ 105 
-57 


NGC4945t 


SBcd 


11? 


u.ui4_ 007 


-20.60 


4.21 


0.30 


+0.45 
-0.18 


NGC5077 


E3 


y 


7 4+ 4 - 7 
'• 4 -3.0 


-25.40 


4.49 


29 


+25 
-13 


NGC5576 


E3 


u 


i-6lg: 3 4 


-24.50 


4.24 


11 


+ 10 

-5 


NGC5813 


E 


y 


6.8l°; 7 


-25.25 


3.97 


20 


+ 17 
-9 


NGC5845 


E3 


n 


2.6±°; 4 


-22.99 


4.25 


2.8 


+2.4 
-1.3 


NGC5846 


E 


y 


nii 


-25.32 


4.21 


24 


+20 
-11 


NGC6086 


E 


y 


°'-ii 


-26.50 


4.42 


78 


+ 07 
-36 


NGC6251 


E2 


y? 


5 0+2.0 

°- y -2.0 


-26.57 


4.73 


96 


+83 
-44 


NGC6264 


S? 


n? 


n 01 +0.004 
u -° -0.004 


-23.61 


3.57 


3.6 


+5.2 
-2.1 


NGC6323 


Sab 


n? 


10+ ' 001 
Uiiu -0.001 


-23.37 


3.36 


2.6 


+3.8 
-1.5 


NGC7052 


E 


y 


7 +2.6 
°-'-1.5 


-25.97 


4.73 


55 


+48 
-26 


NGC7582 


SBab 


u 




-22.83 


3.06 


1.4 


+2.0 

-0.8 


NGC7768t 


E 


y 


13±1 


-26.41 


4.21 


64 


+55 
-30 


UGC3789 


SBab 


n? 


n 1 1+0.005 
u - 1± -0.0O5 


-22.66 


3.18 


1.2 


+ 1.8 

-0.7 



f: Galaxies for which we did not obtain new ARCHANGEL pho- 
tometry. Column (1): Galaxy identifier. Column (2): Morpho- 
logical type. Column (3): Core type, y indicates the galaxy 
contains a core, n indicates the galaxy does not have a depleted 
core. ? indicates that the classification is based on the ve- 
locity dispersion. Column (4): Supermassive black hole mass. 
References are provided in GS13. Column (5): Inclination and 
dust corrected K s absolute spheroid magnitude. For elliptical 
galaxies the typical uncertainty on iC^sph is 0.25 mag. For disk 
galaxies this increases to 0.75 mag, due to the additional dust 
and bulge-to-total corrections. Column (6): (B—K s ) spheroid 
color. Column (7): Spheroid stellar mass. For elliptical galaxies 
the typical uncertainty on M sp h,* is 0.2 dex, for disk galaxies 
this increases to 0.36 dex. 

animation of their surface brightness profiles taken from 
Hubble Space Telescope (HST) imaging. For 19 galaxies 
without HST imaging GS13 assigned a core type based 
on the galaxy's velocity dispersion (see GS13 for further 
details, including a list of those objects with velocity dis- 
persion based core type assignments) . The most massive 
galaxies are typically core-Sersic galaxies and the least 
massive galaxies are exclusively Sersic galaxies, however 
there is a significant region of overlap, with spheroids in 
the 10 10 -10 11 Mq region showing both profile types. For 
all of these galaxies in the overlap region the classification 
was based on their observed light profile. 

2.1. Sources of uncertainty on Mbh and M sp h,* 

The uncertainty on Mbh for each individual object is 
given in Table [TJ these are drawn from the same sources 
as the black hole mass measurements themselves and 
have been adjusted to the distances tabulated in GS13 
- detailed references are given in GS13. The typical 
uncertainty on Mbh is 50%. 

The total uncertainty on M sp h,* has contributions from 
two or three separate components. For elliptical galaxies 
these are the uncertainty on the K s magnitudes mea- 
sured by archangel and the uncertainty in converting 
this magnitude into a stellar mass. The uncertainty in 
the magnitudes derived from the 2MASS photometry are 
determined by the ARCHANGEL pipeline, and are typi- 



cally ~ 0.25 mags (or 0.1 dex). The uncertainty in the 
M/L used to convert these magnitudes to stellar masses 
depends on the uncertainty in the (B—K s ) color, and 
on the uncertain star formation history of each object. 
iBell k de Jong! (|2001h give M/L^ for a range of star 
formation histories, allowing the uncertainty in M/L due 
to the uncertain star formation history to be estimated. 
For our sample, the uncertainty on M/Lr^ (due to both 
the uncertainty on the observed color and the uncertain 
star formation history) is typically 0.17 dex. In ellipti- 
cal galaxies the typical total uncertainty on M sp h,* for 
disk- less systems is 0.2 dex. 

In systems with a stellar disk there is the additional 
source of uncertainty in converting the total magnitude 
into a spheroid magnitude by applying both a dust cor- 
rection and a correction for the bulge-to-disk flux ra- 
tio. The u ncertainty du e to th e dust correction is es- 
timated by iDriver et all (|2008f ) to be 5% in both the 
B- and if -bands. The uncertainty in the bulge-to-disk 
flu x ratio can be est i mate d from the data presented 
in iGraham k Worlevl (|2008l . their table 4), and for the 
galaxies in our sample is typically 0.3 dex. For systems 
with a disk this is the dominant source of uncertainty 
and the total typical uncertainty in M sp h,* is 0.36 dex. 

3. ANALYSIS 

In Figure [3] we show the spheroid stellar mass plotted 
against the supermassive black hole mass for all galax- 
ies in Table Q] We separate galaxies into Sersic (filled 
blue symbols) and core-Sersic (open red symbols). We 
fit separate linear regressions to the Sersic and core-Sersic 
galaxy subsamples, using t he BCES bisector regression of 
lAkritas k Bershadvl ([1996). This technique takes into 
account the measurement uncertainties in both black 
hole mass and stellar spheroid mass and accounts for 
(though does not determine) the intrinsic scatter. For 
the core-Sersic galaxies the best-fitting symmetrical re- 
gression is: 

log ^ = (0. 9 7 ± 0.14) 1 „ g ( J ^ rls -) 

+ (9.27 ±0.09), (3) 

and for Sersic galaxies the best-fitting symmetrical re- 
gression is: 

log = (2.22 ± 0.58) log ( ^ 

Mq ' B \3.Q X 10 10 Mq J 

+ (7.89 ±0.18), (4) 

with rms residuals of 0.47 and 0.90 dex respectively in 
the log Mbh direction. These linear relations are shown 
in Fig. [3]as the solid red and blue lines for the core-Sersic 
and Sersic galaxies respectively. For comparison, the lin- 
ear regression to the combined sample is shown as the 
black dash ed line, which has a slope of 1.50 ± 0.12 (c.f. 
Laor 2001)and an rms residual of 0.67 dex. The best- 
fitting regressions for the two different types of galaxy 
have significantly different slopes that are not consis- 
tent with each other given the confidence intervals on the 
slopes. Sersic galaxies follow an approximately quadratic 
relation, whereas core-Sersic galaxies follow an approx- 
imately linear relation. This is in agreement with the 
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Figure 3. Supermassive black hole mass vs. spheroid stellar mass 
for core-Sersic (open red symbols) and Sersic (filled blue symbols) 
galaxies. The best-fitting linear relations to the two samples are, 
given by Eqns. (J3j and J4{ shown as the solid lines. For compar- 
ison, the best-fitting linear regression for the full sample is shown 
as the dashed line and is dependent on the sample selection. A 
representative error bar is shown in the upper left corner. 



analysis and conclusions of lGrahaml (|2012a|) who studied 
the MBH~M sp h,dyn relation. 

The bend or break in the MBH _ M sp h.* distribution oc- 
curs where the core-Sersic and Sersic relations overlap. 
This is at a spheroid stellar mass M sp h,* ~3x 10 10 M©, 
corresponding to a black hole mass Mbh ~2x 10 s M©. 
The significance of this 'break mass' will be discussed 
in Section [4] Here we simply note that very few Sersic 
galaxies have Mbh greater than this break mass; equally 
no core-Sersic galaxy has a black hole mass below this 
v alue. 

iPark et a l. (2012) examined the effect of using differ- 
ent linear regression techniques on the Mbh _ c relation 
but their results are applicable to supermassive black 
hole scaling relations in general. They reported that 
three popular regression techniques, bces (as used in 
this work) a modified ve rsion of fitexy (jPress et al.l 
ll992tlWemaine etTa l 2002) and a Bayesian technique de- 
veloped bv lKellvl (|2007D . linmix_err return consistent 
results. However, if the measurement uncertainties are 
larger than ~ 15% in the ordinate when using the 'for- 



ward' regression, they find that the BCES routine may be 
biase d to higher slopes (as was noted bv lTremaine et all 
2002). Because of the significant uncertainties on many 
of our low-mass spheroid stellar masses (due to our sta- 
tistical bulge-disk separation) , we redetermined our core- 
Sersic and Sersic relations using the modified fitexy and 
the LlNMix_ERR linear regression methods to test the ro- 
bustness of our result. Following GS13, we determine 
a symmetrical bisector regression from the linmix_err 
routine by determining the 'forward' and 'inverse' regres- 
sions, then determining the line that bisects those two 
regressions. Using the LINMIX_ERR method we find bi- 
sector slopes for the core-Sersic and Sersic samples of 
1.10 ± 0.07 and 2.11 ± 0.46 respectively. We adopt a 
similar approach with fitexy; determining the 'forward' 
and 'inverse' regressions then determining the bisector 
line. This method yields bisector regressions for the core- 
Sersic and Sersic samples with slopes of 1.08 and 2.48 
respectively. With all three symmetric regression meth- 
ods we find consistent best-fitting relations, and in all 
cases the slope for the Serisc galaxies is greater than 2a 
steeper than that for core-Serisc galaxies. We conclude 
that our principal result — that core-Sersic and Sersic 
galaxies follow different MBH _ M sp h,* relations is robust 
against the choice of linear regression method. 

The core-Sersic/Sersic classification is similar to 
the slow rotator/f as t rot ator (SR/FR) classification 
of lEmsellem et "all ( 20071 ) . though as discussed in 
lEmsellem et al.1 (|201lf i" the overlap between the two sys- 
tems is not perfect. Of our objects, 33/75 are part of 
the ATLAS 3D survey (iCappellari et all 120111) and have 
FR/SR classifications from Emscllc m et al.l (2011). We 
also examined the spheroid stellar mass vs. supermassive 
black hole mass relation for these galaxies. We again fit 
linear regressions to the two samples of 9 SR and 24 FR 
galaxies. With this smaller sample we do not find sig- 
nificantly different slopes for the FR and SR samples; 
the slopes for the FR and SR samples are 1.71 ± 0.27 
and 1.32 ± 0.44 respectively. While the slope for FRs is 
somewhat steeper, the difference is not significant given 
the formal uncertainty on the derived slopes. However, 
this may be caused by the galaxies with the most massive 
black holes and the least massive galaxies not having an 
FR/SR classification and therefore not being included 
in these regressions. This same sampling effect, where 
galaxies at the extremes of the supermassive black hole 
mass scaling relations were not well-sampled, led to a 
single log-linear relation being sufficient to describe the 
data in the past. It is only recently, where increased num- 
bers of objects at the extremes of the relations have been 
added, has the bend in the supermassive black hole mass 
scaling relations become evident. Alternatively, the core- 
Sersic/Sersic classification may be more closely linked to 
the mechanism(s) responsible for black hole growth than 
the FR/SR classification. Indeed, FR l enticular galaxies 
with depleted cores are know n to exist (|Dullo fc Grahanl 
l2013tlKrainovic et allf20 13). A larger sample of galaxies 
with both measured supermassive black hole masses and 
kinematic classifications would be desirable. 

From Eqn. (@}, the expected variation of MBH/M sp h,* 
with M sp h,» is close to linear for the Sersic galaxies, 
with the black holes of more massive spheroids repre- 
senting a larger fraction of their host spheroid's mass. 
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Figure 4. The ratio of supermassive black hole mass to spheroid 
stellar mass, Mbh/M sp i,» vs. Mgph *• As previously, core-Sersic 
galaxies are shown as red open symbols and Sersic galaxies as blue 
filled symbols. A representative error bar is shown in the upper 
right corner. 

In Figure H] we show the ratio of supermassive black 
hole mass to spheroid stellar mass as a function of 
spheroid stellar mass. In previous studies which found 



M 



BH 



oc M 



sph,* 



this ratio was thought to be constant, 
with the supermassive black hole ha ving a mass ~ 0.15- 
0.20% of its host spheroid's mass (|Merritt fc Ferraresd 
2001; M arconi fc H unt 200:',;. Here we find that this re- 
mains approximately true for the core-Sersic galaxies, 
albeit with Mbh ~ 0.55% of M sp h,* (the relation for 
core-Sersic galaxies is consistent with a slope of 1 hence 
a constant mass fraction). However, for Sersic galax- 
ies, we find that the average MBH/M sp h,* is offset to a 
lower mean value of 0.3% for our particular sample's mass 
range and it also displays a large range, from ~ 0.02-2%. 

4. DISCUSSION 

4.1. Comparison to previous studies 

In this work we have identified a break in the Mbh~ 
M sp h,* diagram, caused by two separate log-linear re- 
lations with significantly different slopes for core-Sersic 
and Sersic galaxies. While this result is a significant 
change from the commonly accepted view of single log- 
linear scaling relations describing the relationship be- 
tween black holes and their host spheroids, this work is 
no t the firs t to iden tify this change. As noted in Section 
[TJ IGrahaml ([2012a) and GS13 have both previously re- 
ported the need for separate log-linear relations for core- 
Sersic and Sersic galaxies to describe the trend of black 
hole mass with host galaxy dynamical mass and host 
s phero id luminosity respectively. In addition, IGrahaml 



sp. 
(2( 



Mbh _ L relations cannot both be log-linear due to the 
non-linear L-cr relation for early-type galaxies. 

While in this paper we have quantified the first bent re- 
lation between black hole mass and host spheroid stellar 
mass, we can compare our results with past work per- 
taining to the host's dynamical mass. GS13 reported an 
Mbh-L^s relation for the Sersic galaxies with a slope of 
2.73 ± 0.55, which they equated to a slope of 2.34 ± 0.47 



in the MsH-Mdyn diagram using M/Lk, 



oc 
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2007) had previously noted that both the Mbh-c and 



(e.g., 

lLa Barbera et al.ll2010t iMagoulas et al.ll2012| ). Prior to 
that, Graham (2012) had reported a slope of 1.92 ± 0.38 
for the MBH-Mdyn relation based on independent data. 
These two steep slopes compare well with our slope of 
2.22 ± 0.58 for the Sersic galaxy MBH-M sp h,* relation. 
For the core-Sersic galaxies, GS13 had reported a slope 
of 1.10 ± 0.20 in the Mbh _ La" s diagram. If core-Sersic 
galaxies are predominantly built by dry mergers of galax- 
ies near or above the high- mass end of the Sersic distribu- 
tion, then the slope in the MBH~M sp h.* diagram should 
be the same as the slope in the Mbh _ L diagram. Graham 
(2012) had additionally reported a slope of 1.01 ± 0.52 
for the core-Sersic galaxies in the MBH _ Md yn ,sph diagram 
(with the large uncertainty reflecting their small sample 
size) . These two shallower slopes are consistent with our 
slope of 0.97 ± 0.14 in the MBH~M sp h,* diagram. 

Other authors have noted that low luminosity (or 
mass) galaxies are consistently offset from single log- 
linear black hole scaling relations derived from samples 
dominated by massive systems. iGreene et al.l (2008) 
identified a population of low mass galaxies (stellar 
masses in the range 10 9 -10 10 M Q ) whose black hole 
masses were o ffset below the MBH _ Mb u igc relation of 
lHaring fc Rbd (I2004T) by an or der of magnitude. More 
recently. IMathur et al.l (|2012D determined black hole 
masses for a sample of 10 narrow-line Seyfert galaxies 
with host spheroid luminosities in the range 3 x 10 — 
3 x 10 10 L Q (corresponding to a range in stellar mass 
of approximately 10 10 -10 11 Mp,) and again found their 
galaxies to be offset below the lHaring fc Rixi (|2004D rela- 
tion. While IMathur et al.l ([20121 attribute much of this 
offset to 5 of their objects being 'pseudobulges', the re- 
maining 5 are 'classical' bulges and are st ill substantially 
offset fr om the Mbh - Lhnip . relation of iGiiltekin et all 
(|2009fl . IGreene et all ((2008) also identify many of their 
objects which are significantly offset below the classical 
log-linear MsH^Mbulge relation as being well-fit with a 
de Vaucouleurs (n = 4) profile, suggesting they are not 
pseudobulges. The offset objects found by both these 
studies are consistent with the log-linea r scaling relations 
for Sersic galaxies reported in this work, IGrahaml (|2012af) 
and in GS13. While we do not identify pseudobulges in 
our sample, we note that both NGC4486a and NGC821, 
both 'classical' elliptical galaxies with no indication of 
a pseudobulge, are consistent with our Sersic scaling re- 
lation and offset from the core-Sersic relation, arguing 
against the offset nature of the low-mass systems being 
a pseudobulge phen omena. 

iSani et al.l poll . iVika et al.l (pOll and lBeifiori eTall 
(l20ll " have all recently constructed MgH -Msph,* rela- 
tions (or in the case of IVika et all 120121 the closely- 
related MBH-Lif,s P h relation) for large samples (~ 50 
galaxies). All three studies only considered single log- 
linear fits to their data and are dominated by massive 
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galaxies with M sp h.* > 3 x 10 10 M©. As expected from 
their high-mass-dominated samples, all three studies find 
MsH-Msph,* relations with slopes ~ 1, consistent with 
our finding for core-Sersic galaxies. However, in all three 
studies a number of galaxies are offset below those au- 
thors' single log-linear relations. With host spheroid 
masses around 3 x 10 10 M©, they are consistent with the 
high-mass regime of our Sersic galaxy relation. These 
studies highlight the need to extend the range of super- 
massive black hole masses and host spheroid masses used 
to examine the black hole scaling relations. 

4.2. Scatter about the MBH - Miph,* relation 

As well as determining the slope of the supermassive 
black hole scaling relations, many studies also examine 
the scatter about the relations. In particular, this is done 
to argue that one of the relations has reduced intrinsic 
scatter compared to the other common scaling relations, 
and therefore is the 'fundamental' driving relation. At 
low masses the scatter in our Figure [3] is dominated by 
the large measurement uncertainties due to our statisti- 
cal dust and bulge correction, making such a comparison 
difficult for this study. Instead, we will make a few simple 
observations on the revised expectations for the intrinsic 
scatter as a result of our bent scaling relation. 

As noted in Section (3] the observed scatter (the sum 
of the intrinsic scatter and measurement uncertainty) in 
the black hole mass direction is significantly larger for the 
Sersic galaxies than for the core-Sersic galaxies: 0.90 dex 
compared to 0.47 dex. An increase is expected given the 
typically larger measurement uncertainties for the Sersic 
galaxies due to the statistical dust and bulge correction. 
However, we also expect the scatter in the vertical direc- 
tion to increase because of the increased slope of the re- 
lation for the Sersic galaxies. This is not the case for the 
intrinsic scatter orthogonal to the relation — we would 
expect this to be reduced relative to a single log-linear 
relation fit to the entire data. The orthogonal scatter 
we find for the core-Sersic galaxies is 0.26 dex, and for 
the Serisc galaxies it is 0.29 dex, an improvement over 
the 0.31 dex orthogonal scatter we find for the single 
log-linear relation. 

Finally, a number of theoretical studies have examined 
the idea of the supermassive black hole scaling relations 
being the product of the repeated merging of a randomly 
see ded initial population of black holes and host galax- 
ies (|Pengl 120071 iJahnke fc Macciol 120111 ). These studies 
predict a decrease in t he intrinsic scatter as h ost galaxy 
mass increases (though Jahnke & Maccio 2011, note that 
the addition of star formation to their model reduces the 
rate of this decrease in scatter) . Given the large measure- 
ment uncertainties for many of our Sersic galaxies it is 
difficult to quantify any variation of the intrinsic scatter 
with host mass in our sample. Moreover, we expect that 
only for the core-Sersic galaxies is the binary merging of 
both galaxies and black holes the dominant mechanism of 
growth. As these galaxies only span a relatively narrow 
range in host spheroid mass, any systematic variation of 
the intrinsic scatter within this subsample is uncertain. 

4.3. Comparing supermassive black hole and nuclear 
star cluster scaling relations 

iFerrarese etaLl (|2006bD and iWehner fc Harris! (|2006fi 
have argued that nuclear star clusters and super- 




Spheroid stellar mass [M ] 

Figure 5. Central massive object mass, Mcmo vs - spheroid stel- 
lar mass, Mgpi^* for our Sersic galaxies (blue points) and for the nu- 
clear star clusters (and host spheroids) of Scott & Graham] 1 12013 , 
open green squares). The blue and green lines are the best-fitting 
symmetric al linear relations to the two samples (the green line is 
taken from lScott fc Graham! i2013ft ). 

massive black holes form a single class of cen- 
tral massive object (CMO) based on their al- 
legedly c ommon mass scalin g relations, though recent 
studies (IBalcells et all 120071: IGraham fc Spitlerl 1 2009 1 : 
| Grahaidl2012bHErwin fc Gadottill2012l:lLeigh et alj20ia 
iScott fc Graham] I2013D have since argued against this 
scenario. We briefly re-examine the connection between 
nuclear star cluster and supermassive black hole scaling 
relations in the light of the bent MBH-M sp h,* relation. In 
Figure[S]we show Mcmo vs. M sp h,* for the Sersic galax- 
ies presente d in this study and the nuclear star clusters 
presented in lScott fc Graham! (]2013l ). The two lines show 
the best-fitting linear relations to the two samples. The 
MNC - M S ph,* line is taken from IScott fc Graham! ()2013f ) 
and has a slope of 0.88±0.19. The relation for the nuclear 
star clusters is 2.3a shallower than that for the super- 
massive black holes. T his difference is m o re pro nounced 
than that reported by IScott fc Graham] (|2013l ) , due to 
our improved bent supermassive black hole scaling rela- 
tion, strongly suggesting that supermassive black holes 
and nuclear star clusters do not follow a common scaling 
relation and therefore do not have a common formation 
mechanism. 

4.4. The growth of supermassive black holes and stellar 

spheroids 

As noted in Section [3] the break or bend in the 
MBH-M S ph,* relation occurs at M sp h,* — 3 x 10 10 M©. 
This mass scale has been identified by a number of au- 
thors as marking a change in early-type galaxy proper- 
ties. Below this mass g alaxies are typically you ng, low- 
surface density objects (Kauffm ann et al.ll2003al ). whose 
spheroid's surface brightness is well-fit by a Sersic pro- 
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file (|Graham fc Guzmanl [200l iBalcells et al.l 12003ft . and 
define a sequence of increasing M/L and bulge fraction 
(| C app ellaril 1 2 1 2ft — until dwarf spheroidal systems ap- 
pear around Mb ^ — 14 mag. Above this mass galaxies 
are typically old with high Sers ic indices (jGraham et al.1 
12001c iKauffmann et aD l2003afl , have core-Sersic sur- 
face brightness profiles (|Graham fc Guzmanl I2003D . are 
bulge-dominated an d show only a narrow r ange in M/L 
(jCappellaril 120121) . iTremonti et~all (|2004l) also report 
that the mass-metallicity relation flattens above this 
mass scale. 

The results discussed above are all consistent with 
a scenario in which galaxies with spheroid masses be- 
low the rough 3 x 10 10 M Q divide grow predominantly 
through gas-rich processes involving significant in-situ 
star formation whereas the growth of those galaxies 
above this mass scale is dominated by 'dry' processes 
involving little gas and the direct accretion of stars 
via mergers. This i s similar to a scenario adopted by 
IShankar et al.l ((201 2f ) based on their semi-analytic model 
of black hole and galaxy formation. 

In the high-mass regime where little gas is involved, 
black hole growth will be dominated by the binary 
merger of supermassive black holes during galaxy merg- 
ers. In a dry merger, to first order, the masses of the 
two progenitor galaxies will simply be added, as will 
the masses of their supermassive black holes. Thus, in 
this regime we can roughly expect a linear relationship 
between supermassive black hole mass and host stellar 
mass, as observed. 

In the low mass regime gas plays a more significant 
role, complicating the picture. Galaxies can grow from 
the direct cooling of gas from their immediate surround- 
ings, and even when mergers do play a significant role 
the resulting growth (of both stellar mass and black hole 
mass) will not be simply additive as the 'wet' mergers will 
likely be accompanied by significant star formation. The 
slope of ~ 2.2 in the MBH-M sp h,* relation that we find for 
Sersic galaxies implies that black holes have grown more 
rapidly than their host spheroids. This is consistent with 
the growth of these galaxies being dominated by gaseous 
processes (e.g., gas-ri ch mergers, cooling of hot ga s, sec- 
ular evolution: see Hopkins fc Hernquistl |2009, for a 
detailed discussion of the full range of processes) that 
efficiently channel material on to the central black hole 
(see e.g., iMarconi et al.l 12004ft but are less effective at 
growing the host galaxy spheroid. One such example of 
a supermassive black hole growing more rapidly than its 
host g alaxy has recently been observed bv lSevmour et ahl 
(2012). This steepening of the supermassive black hole 
mass - host spheroid stellar mass r elation has been repro- 
duced in a number of simulations (|Cirasuolo et al.ll2005l : 
iDubois et afll20Tl iNeistein fc Netzerll2013l figures 5, 3 
and 8 respectively), though these authors have typically 
focused on the agreement with the classical single log- 
linear relation in massive galaxies, rather than the down- 
turn at lower mass in their simulated data. 

5. CONCLUSIONS 

The results in Section [3] represent a major revision to 
the accepted picture of the relationship between super- 
massive black holes and their host galaxies. Our results 
have significant implications for: studies of how super- 
massive black holes and their hosts co-evolve over cosmic 



time; the detection of intermediate-mass supermassive 
black holes; and models of feedback from active galactic 
nuclei due to supermassive black hole growth. 

We have examined the bent nature of the supermas- 
sive black hole mass - spheroid stellar mass relation for 
a large sample of 75 galaxies. We find that, when sepa- 
rating galaxies into core-Sersic and Sersic galaxies based 
on their central light profiles, they follow significantly 
different supermassive black hole scaling relations whose 
slopes are not consistent with one another. For the core- 
Sersic galaxies Mbh oc M° p ^ 7 ^ 014 , whereas for the Sersic 
galaxies Mbh oc M 2 ^ 2 ^ ' 58 . These results are consistent 
with the expectation from a single log-linear supermas- 
sive black hole mass - host velocity dispersion relation 
combined with a bent relationship between host spheroid 
luminosity and host velocity dispersion. We find that 
the supermassive black hole mass is an approximately 
constant 0.55% of a core-Sersic galaxy's spheroid stellar 
mass. The non-log-linear MBH~M sp h,* relation implies 
that, for Sersic galaxies, Mbh is not a constant fraction 
of the host spheroid mass, but represents an increasing 
fraction with increasing M sp h.*. 
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